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^ ' Abstract 

A scale-invariant spectrum of isocurvature perturbations is generated during col- 
lapse in the ekpyrotic scaling solution in models where multiple fields have steep neg- 
i— i" ative exponential potentials. The scale invariance of the spectrum is realized by a 

■ tachyonic instability in the isocurvature field. This instability drives the scaling solu- 
tion to the late time attractor that is the old ekpyrotic collapse dominated by a single 

' field. We show that the transition from the scaling solution to the single field domi- 

nated ekpyrotic collapse automatically converts the initial isocurvature perturbations 
about the scaling solution to comoving curvature perturbations about the late-time 
attractor. The final amplitude of the comoving curvature perturbation is determined 

■ by the Hubble scale at the transition. 

in ! 

1 Introduction 

[ — | The existence of an almost scale-invariant spectrum of primordial curvature perturbations on 
large scales is one of the most important observations that any model for the early universe 
should explain. An inflationary expansion in the very early universe is most commonly 
^ . assumed to achieve this, but it is important to consider whether there is any alternative 
model. In this paper we focus on the ekpyrotic scenario as an alternative [1] (see also [2j [3]). 
In the old ekpyrotic scenario, the large scale perturbations are supposed to be generated 
during a collapse driven by a single scalar field with a steep negative exponential potential. 
It was shown that the Newtonian potential acquires a scale-invariant spectrum, but the 
comoving curvature perturbation has a steep blue spectrum [I]. In this scenario we need a 
mechanism to convert contraction to expansion, and for a regular four-dimensional bounce, 
the scale-invariant Newtonian potential is matched to the decaying mode in an expanding 
universe, and the growing mode of curvature perturbations acquires a steep blue spectrum 
[S, E] . It has been suggested that this conclusion might be altered by allowing a singular 
matching between collapse and expanion [7j, but the general rule that the comoving curvature 
perturbation remains constant still holds for adiabatic perturbations on large scales [8]. In 
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a braneworld context, the conversion from contraction to expansion might be accomplished 
by a collision of two branes where one of extra-dimensions disappears [9]. It was argued 
that the scale-invariant Newtonian potential can be transferred to the comoving curvature 
perturbations by this singular bounce |10| . However, without having a concrete theory to 
describe the singularity, it is difficult to have a definite conclusion on how perturbations pass 
through the singularity. 

Recently, there has been some progress in generating a scale invariant spectrum for 
curvature perturbations in the ekpyrotic scenario [TTI [T2"| IT3"]. by considering non-adiabatic 
perturbations which has been suggested previously by Ref. [15]. In this case we require two 
or more fields. If these have steep exponential potentials then there exists a scaling solution 
where the energy density of the fields grow at the same rate during collapse [HH [T7] . The 
isocurvature perturbations then have a scale-invariant spectrum [IB]. These isocurvature 
perturbations can be converted to curvature perturbations if there is a sharp turn in the 
trajectory in field space [TTJ [121 [13] . F° r example Ref. [11] considered a situation where one 
of the fields changes its direction in field space, which corresponds to a time when a negative 
tension brane is reflected by a curvature singularity in the bulk, in the context of the heterotic 
M-theory. Refs. [121 E3] considered a regular bounce realized by a ghost condensate. One 
of the fields exits the ekpyrotic phase and hits the transition to the ghost condensate phase 
that creates a sharp turn in the trajectory in field space and curvature perturbations can be 
generated [12]. It is still necessary to match curvature perturbations in a contracting phase 
to those in an expanding universe, but it is shown that the comoving curvature perturbation 
is conserved on large scales resulting in an almost scale-invariant spectrum observed today 
for a regular bounce like a ghost condensate model [T2] . 

The isocurvature perturbations behave like Ss oc H on large scales. As the Hubble 
parameter is rapidly increasing in a collapsing universe, this signals an instability. In fact it 
is easy to see that we always require an instability of this form in order to generate a scale- 
invariant spectrum with a canonical scalar field during collapse^. As the amplitude of field 
perturbations at Hubble exit are of order H, we require the super-Hubble perturbations to 
grow at the same rate to maintain a scale-invariant spectrum. Ref. [TJJ] studied this instability 
in a phase space analysis. The multi-field scaling solution was shown to be a saddle point in 
field space and the late-time attractor is the old ekpyrotic collapse dominated by a single- 
field. A tachyonic instability drives the scaling solution towards the late-time attractor (see 
also HH[I2]). 

In Ref. [12], we pointed out that the natural turning point in the field space trajectory 
due to the instability of the scaling solution might itself offer the possibility of converting the 
scale invariant spectrum of isocurvature field perturbations into a scale invariant spectrum 
of curvature perturbations. In this paper, we confirm this expectation by explicitly solving 
the evolution equations for perturbations in a two field model. We find that the ratio 
between curvature perturbations and isocurvature perturbations at the final old ekpyrotic 
phase is solely determined by the ratio of exponents of the two exponential potentials, and 
the amplitude is set by the Hubble rate at the transition time. 

4 The only way to produce a scale-invariant spectrum without the presence of an instability seems to 
be due to non-canonical kinetic terms, as in the case of axion-type fields which can acquire scale-invariant 
perturbation spectra while remaining massless [18) . 
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2 Homogeneous field dynamics 



We first review the background dynamics of the fields. During the ekpyrotic collapse the 
contraction of the universe is assumed to be described by a 4D Friedmann equation in the 
Einstein frame with scalar fields with negative exponential potentials 

3H 2 = V + ^l (1) 

where 

V = -YtVie-*** , (2) 

i 

and we take Vi > and set 8nG equal to unity. 

The authors of [H] found a scaling solution (previously studied in [T6], HZ]) in which both 
fields roll down their potential as the universe approaches a big crunch singularity. In this 
ekpyrotic scaling collapse we find a power-law solution for the scale factor 

2 1 

a oc (-t) p , where p = ~ < ~ > ( 3 ) 



where 



2 ~~ -Vje-^i ~ cf 



(4) 



As we will see in the next section, it is possible to generate scale-invariant isocurvature per- 
turbations around this background. However, the ekpyrotic scaling solution PJ is unstable. 

In addition to the scaling solution we have fixed points corresponding to any one of the 
original fields (pi dominating the energy density where the other fields have negligible energy 
density. These correspond to the original ekpyrotic power-law solutions where 

2 

a oc (— t) Pi , where Pi = , (5) 

for cf > 6. We find that any of these single field dominated solutions is a stable local 
attractor at late times during collapse. 

In Ref. [12], the ekpyrotic scaling solution (J3J) was shown to be a saddle point in the 
phase space. We briefly review the phase space analysis. Introducing phase space variables 

psaisnin] 

the first order evolution equations for the phase space variables are given by 

^ = -^{l-^-c^-yl (8) 

jfi = Vi ( 3 E^- C J^) > ( 9 ) 
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where N = log a. The Friedmann equation gives a constraint 

3 j 

There are (n + 2) fixed points of the system where dxi/dN = dyi/dN = 0. 



A : 



0. 



(10) 
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B 
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3p c. 



= 0, (for j ^ z) 



A 




— i 



(12) 



(13) 



In this paper, we focus on the fixed points B and Bi assuming c\ > 6 and c ~i 2 < 1/6. The 
linearized analysis shows that the multi-field scaling solution, B, always has one unstable 
mode. On the other hand, the single field dominated fixed points, Bi, are always stable. 

From now on we concentrate our attention on two fields case. Then we have three fixed 
points B, B\ and B 2 . It is interesting to note that in the (x±, x 2 ) plane, the fixed points B, 
B\ and B 2 are connected by a straight line, which is given by 

ClC 2 



C 2 Xi + C\X 2 



(14) 



The eigenvector associated with the unstable mode around the scaling solution B lies in the 
same direction as the line (fT4l) . Thus this is an attractor trajectory, which all solutions near 
B approach. Fig. 1 shows numerical solutions for the evolution of x\. Initial positions in the 
phase space are perturbed away from B along the line (Tl4l . The solutions go to B\ or B 2 , 
depending on the initial position in the phase space. 

An important observation is that, as we follow phase space trajectories during the tran- 
sition from the scaling solution B to the attractor solutions Bi or B 2 , the solutions obey the 
relation Eq. (|14|) . even far away from the saddle point, B. Using the Friedmann equation 
and the field equations, we can show that 



H 



Cl 








and hence 



h - ti _ h. = — 

Ci c 2 a 3 ' 



(15) 



(16) 



where C is an integration constant. In terms of <p\ and (p 2 , Eq. ( IT41) can be rewritten as 

^ + ^ = H. (17) 

Cl c 2 

and hence we see that for trajectories starting from point B we have C — 0, which is the 
late time attractor. Thus we see that Eq. fflTj) holds even during the transition caused by 
the tachyonic instability from point B to B\ or B 2 and in the final single-field dominated 
phase. This fact will be important when we study perturbations. 

We will study the behaviour of perturbations during the transition in the next section. 
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Figure 1: Left: Numerical solutions for x±(N). The horizontal axis is iV = log a and we take 
Ci = 40 and C2 = 30. The initial time is N = 0.05. Note that N decreases towards the future 
in a collapsing universe. Right: The corresponding phase space trajectories in (xi,X2) plane. 



3 Generation of quantum fluctuations 



In this section, we consider inhomogeneous linear perturbations around the background 
solution. We consider the scalar field perturbations on spatially flat hypersurfaces. Then 
the scalar field perturbations are given by [22j [231 E 



k 2 

+ + —I 



c 2 V i exp(-c i (j) i )6(j) i - ^ ( jj<fii< 



0. 



(18) 



We can decompose the perturbations into the instantaneous adiabatic and entropy field 
perturbations as follows [24] : 



Sr 



+ 



6s 



(19) 



The adiabatic field perturbation Sr is the component of the two-field perturbation along the 
direction of the background fields' evolution while the entropy perturbation 6s represents 
fluctuations orthogonal to the background classical trajectory. The adiabatic field perturba- 
tion leads to a perturbation in the comoving curavture perturbation: 



Tic 



H6r 



(20) 



whereas the entropy field perturbations correspond to isocurvature perturbations. 
Their evolution equations are given by [24] 



k 2 

6r + 3H6r H — -Sr + 
a 2 



Vr 



e 2 - 



a 3 r 2 " 



H 



6r = 266s + 2 



MM' 



6s, (21) 



k 2 

5s + 3H5s + —5s + (V ss - 9 2 )5s 



r5r 



IE 



+ r 5r 



where the angle 9 is defined as 

cos 9 

such that 



sin 9 



r 
9 



(cos 
V, 



2 + (sin 



l) 



and 



Vs., 



(sin^)ciVi exp(-ci0i) + (cos6 , )c 2 l / 2exp(-C202), 
(cos^)ciVi exp(-ci0i) - (sin6 ) )c2V 2 exp(-c 2 02), 

— (sin 9) 2 clVi exp(— ci<pi) — (cos 9) 2 c^V2 exp(— c 2 4>2) 

— (sin 9) 2 c^V2 exp(— c 2 2 ) — (cos 9) 2 c 2 Vi exp(— Ci0i) 



For the multi-field scaling solution, B, we have 

9 = arctan — 



and for the single field scaling solutions we have 



9 



7T 



- (B 1 ), 9 = (B 2 ). 



(22) 

(23) 

(24) 
(25) 



(26) 
(27) 
(28) 
(29) 



(30) 



(31) 



Thus we have 9 =constant for the fixed points and the adiabatic and the entropy fields are 
decoupled. This allows us to quantise the independent fluctuations in the two fields. 

For the multi-field scaling solution B, the spectrum of quantum fluctuations of the entropy 
field is given on large scales (k <C aH) by 



V, 



_fc3 
2^T 



.\5s 2 \ = Cl k - 2 {-krf-'" 
a 1 



where r < is conformal time, and 
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3e 



H/H 2 = l/p, 



4 (e-1, 

and C v = 2 u ^ 2 T(u)/tt 3/2 [EJ. The spectral tilt is given by 

2 

An*, ~ - , 



(32) 



(33) 



(34) 



to leading order in a fast-roll expansion (e ^> 1) [TH [121 [13] . in this limit, the spectrum ([321 
can be written as 

Vll 2 = e ^ . (35) 



2tt 



Note that \H\ is rapidly increasing and thus Ss is also growing on super-Hubble scales due to 
the tachyonic instability This instability is essential in order to realize the scale invariance 
of the spectrum (I35I) . The amplitude of field perturbations at Hubble exit is of order H 
and thus we require the super-Hubble perturbations to grow at the same rate in order to 
maintain a scale-invariant spectrum. 

Note that in the simplest model, the spectrum is slightly blue [HI [121 13] • However, any 
deviations from an exponential potential for adiabatic field can introduce the corrections to 
the spectral tilt and thus it becomes model dependent [ITJ IT2] . 

The spectrum of quantum fluctuations in the adiabatic field about the scaling solution 
has the same power-law form on large scales 

V Sr = C^(-kr) 1 - 2 ^ (36) 

p or 

where to fi ~ 1/2 to leading order in 1/e. Thus the adiabatic field perturbations become 
constant in the large scale limit and the spectral tilt is given by 

An Sr ~ 2 . (37) 

Thus we have 

^ oc {-krf , (38) 

and hence in what follows we can neglect the adiabatic field fluctuations in the large-scale 
limit. 

By contrast, for the single field dominated scaling solutions, the adiabatic and entropy 
field perturbations are both frozen on super-Hubble scales: 

Ss , 5r — const. (39) 

These perturbations both have a steep blue spectrum if they cross the horizon when the 
background solutions are described by the single field dominated solution [T9~] . 



4 Generation of curvature perturbations 

Now let us consider the evolution of perturbations in a situation where the classical solution 
starts from near the saddle point, B. As emphasized in Ref. [19J this requires an additional 
preceding mechanism that drives the classical background solution to the unstable saddle 
point throughout our observable part of the universe. In this paper, we will not discuss the 
mechanism required to bring the classical solution to the saddle point and we just assume 
that the classical solution stays near the saddle point for long enough to ensure that a scale- 
invariant spectrum of isocurvature perturbations is generated over the relevant scales for the 
observed large scale structure of our Universe. 

Then the initial conditions for the adiabatic and the entropy field perturbations can be 
set from the amplitude of quantum fluctuation as described in the previous section: 

5r = 0, Ss = e — , (40) 
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on sufficiently large scales and for 1/e < 1. 

Unless the spatially homogeneous background solution is located exactly at the fixed 
point, the tachyonic instability drives the background solution away from the multi-field 
scaling solution, B, to one of the single field dominated solutions, B\ or B 2 , depending on 
the initial conditions. During the transition 8 is not constant and the adiabatic and entropy 
field perturbations mix, so it is possible to generate perturbations in the adiabatic field, and 
hence comoving curvature perturbations (I20p . from initial fluctuations in the entropy field. 

We can solve the evolution equations (12 ip and (1221) numerically for any given classical 
background. Figure 2 shows the behaviour of Sr and 5s. Due to the coupling between Sr 
and 5s during the transition, curvature perturbations are generated during the transition. 
On the other hand, 5s shows a tachyonic instability according to Eq. fHUl) close to the 
scaling solution, but when the background solution goes to B iy the entropy field perturbation 
becomes constant. The final amplitude of 5r depends on when the transition from B to Bi 
occurs, but, interestingly, the final ratio between 5r and 5s does not depend on the details 
of the transition. We find that the ratio is determined solely by the parameters c\ and c 2 as 

<)r " at Ex, (41) 



5s 
Sr 
5s~ 



Cl 
C 2 



C2 



— , at Bo 



(42) 



as is shown in Fig. 2. We will explain later why such a simple result is found. 
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Figure 2: Left: Solutions for 5r(N), using the same parameters as in Figure 1. The corre- 
sponding background solutions are shown in Figure 1. Right: The ratio between 5r and 5s. 
The ratio approaches a constant given by Eqs. fj4Tj) and (142]) . In this case 1.3333 for B\ and 
-0.75 for B 2 . 

The resulting curvature perturbation on a comoving hypersurface in the final single-field 
dominated phase is thus given by 

at Bi, (43) 
at B 2 . (44) 



—Sr - — Sr. 

Cl 



H 



-Sr = — Sr, 

i 2 °2 
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The equation for the entropy field perturbation ( |22i) can be rewritten as [21] 



5s + 3H5s + ( % + V„ + 36 2 ) 5s = 4-^* , (45) 
\ a 1 j r cr 

where ^ is the curvature perturbations in Newtonian gauge. The change of the curvature 
perturbations is determined by [24] 

U c = S-^ + — 65s. (46) 
Ha r 

On large scales, we can neglect (k 2 /a 2 )^ and it is possible to reproduce the previous results 
from Eq. (I46| . 

Although the physical meaning of the instantaneous adiabatic and entropy field perturba- 
tions is clear, the dynamics of the perturbations during the transition are rather complicated 
in this basis. We find it is much easier to work in terms of new variables [19] 



_ C 2 01 + Ci0 2 ^ _ Ci0i - C 2 02 

\jc\ + c\ \fc\ + c\ 

corresponding to a fixed rotation in field space. The potential Eq. ([2]) can then be simply 
re-written as [291 [TBI [19] 

V = -U( X )e~ ap , (48) 

where 

\ = Y.\> (49) 
and the potential for the orthogonal field is given by 

- U{ X ) = -V x e^ Cl/c2)cx - V 2 e (c2/ci)cx , (50) 

which has a maximum at 



cVA 

c 2 + c 2-\c 2 2 V 2/ 

The multi-field scaling solution corresponds to x = Xoj while ip is rolling down the 
exponential potential. The potential for x has a negative mass-squared around x = Xo, an d 
thus x represents the instability direction. If the initial condition for x is slightly different 
from xo or X is not zero, then x starts rolling down the potential and the solution approaches 
a single-field dominated solution. 

Note that perturbations 5(p and 5x coincide with the instantaeous adiabatic and entropy 
field perturbations respectively, defined in Eq. f|T9|) . at the scaling solution, B. Thus we 
can use the initial perturbations (I4"01 due to vacuum fluctuations about the scaling solution 
previously calculated. However as we follow the evolution away from this saddle point we 
can no longer identify tp and x with the adiabatic and entropy perturbations. Nevertheless, 
the dynamics of perturbations turns out to be much simpler using these fields. 
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In terms of tp and X the equations for perturbations are given by 

k 2 

5^ + 3H5(p + —dip + M^Sip + M VX 5 X = 0, (52) 
k 2 

5 X + 3H5 X + — S X + M XX 5 X + M^Sip = 0, (53) 

d 



where 



M w = V w - 1 ^0 2 ) , (54) 
M *x = V ,m ~ ^ (^x) , (55) 

^xx = ^xx - ^ (^X 2 ) • (56) 

A key observation is that the phase space trajectory of the background fields during the 
transition, Eq. (fl7|) . can be re- written as 

| = * (57) 

Thus even away from the multi-field scaling solution, ip obeys a scaling relation. We can 
then show that two of the effective mass terms become 

M w = V }VV + cV t y = 0, (58) 
M vx = V vx + cV x = 0, (59) 

independently of the form of U(x) since V oc exp(— ap). 

Thus on large scales 5ip is constant. As we take Sip = as our initial condition, this 
remains so even during the transition and in the final single field dominated phase. Thus 
from Eq. (|47p we have a relation between <50i and <50 2 , 

c 2 <50i + Cl 5<p 2 = 0, (60) 

and and (50 2 are given by 

n S X , (61) 



r 2 
c 2 



S X . (62) 



c? + c| 



In the single field dominated solutions, f is no longer the adiabatic field. The adiabatic field 
is simply 0i (or 2 ) in B\ (or S 2 ) and the entropy field is — 2 (or 0i) in I?i (or B 2 ). Thus 
at the late-time attractor we have 

Sr = 5(f>i, 6 s = — <50 2 , at B\, (63) 
<5r = 502, 5s = (50i, at £> 2 . (64) 
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Thus the ratio between 5r and 5s is determined by the ratio between 8<pi and 5<p2 and we 
can easily find the ratio 5r/5s as given in Eqs. fHTj) and (H2"l) . 

The amplitude of the curvature perturbations can be estimated from 5\- In the initial 
stage, 5x coincides with the entropy field perturbations 5s, and thus its initial amplitude on 
super-Hubble scales is given by 



H 



2tt 



(65) 



where we used e = c 2 /2 in Eq. (HOl) and c 2 is given by Eq. (H91 . After the transition, 5\ 
becomes massless and the amplitude becomes frozen on large scales. In the single field 
dominated solutions, the comoving curvature perturbation 1Z C is given by 



1 



c\ + cl 



(66) 



Assuming the transition occurs suddenly, the final amplitude of the comoving curvature 
perturbation is given by 



cf + c| 



H 



2vr 



(67) 



where the subscript T means that the quantity is evaluated at the transition time. On the 
other hand the amplitude of the entropy field perturbation is given by 



5s 



5s 



c\ + cl 



H 



cl + c 2 2 



2tt 
H 



2tt 



at B%, 



at Bo 



T 



(68) 
(69) 



From the numerical solutions for 5r and 5s, we can reconstruct Ht- We confirm that Ht 
constructed in this way agrees with the Hubble scale at the transition as is shown in Figure 3. 



5 Conclusion 

In this paper we have studied the generation of curvature perturbations during an ekpyrotic 
collapse with multiple fields. We must assume that the classical background solution starts 
from a state very close to a saddle point in the phase space that corresponds to the multi- 
field ekpyrotic scaling solution. If the solution stays at this fixed point for long enough, a 
scale-invariant spectrum of isocurvature perturbations is generated over the range of scales 
that is relevant for large scale structure in the Universe. Even if the background solution 
deviates only slightly from the multi-field scaling solution initially, the tachyonic instability 
eventually drives the solution to the old ekpyrotic collapse dominated by a single field. 
During this transition, the initial isocurvature field perturbations generate a scale-invariant 
spectrum of comoving curvature perturbations. 

First we studied the perturbations by decomposing them into the instantaneous adiabatic 
and the entropy field perturbations. These fields are decoupled at the fixed points. The adia- 
batic field perturbations are effectively massless around both the multi-field scaling solutions 
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Figure 3: Left: Solutions for Sx{N). We used the same parameters as Figures 1 and 2. Right: 
Solutions for log \ H\ with three different background solutions. We also show log \ Ht\ that 
is determined from the numerical solutions for 8\- 



and the single field dominated solution, so become constant on large scales. On the other 
hand, the entropy field perturbations have a tachyonic mass around the multi-field scaling 
solution, where they grow like 5s oc H on large scales, but they are effectively massless 
around the single field dominated solution. We set initial conditions for these perturbations 
from quantum fluctuations about the multi-field scaling solution. As the adiabatic field 
perturbations have a blue spectrum they can be neglected compared with the entropy field 
perturbations on large scales. However, during the transition to the single-field attractor, 
adiabatic field perturbations are generated from the entropy field perturbations. 

It turns out to be more convenient to use new fields ip and x defined in Eqs. (1471) to 
follow the perturbations through the transition. These fields coincide with the adiabatic and 
entropy fields around the multi-field scaling solution but not during the transition or at the 
final single field dominated phase. In terms of <p and Xi t ne potential is given by a product of 
the potential for x, U(x), and an exponential potential for (p. U(x) has an extremum which 
corresponds to the multi-field scaling solution. A crucial observation is that even during the 
transition and in the final phase, ip satisfies the scaling relation Eq. (157)) . We can then show 
that the field perturbations for <p> and x are always decoupled. Since 5ip = on large scales 
during the ekpyrotic scaling solution, this remains so. This determines the ratio between the 
adiabatic and entropy field perturbations at the final phase. On the other hand, 5x grows 
during the scaling solution and then becomes constant during the single field dominated 
solution. 

Our final results are Eqs. (ISTj) . (168)) and (169)) for the amplitude of comoving curvature 
and isocurvature field perturbations during the single field ekpyrotic collapse phase. The 
amplitude of the comoving curvature perturbation is determined by the Hubble scale at the 
transition. 
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We still need to convert the ekpyrotic collapse to expansion (see, for instance, [T2l [T3] ) 
and see how this curvature perturbation is matched to that in an expanding universe. For 
a regular bounce the comoving curvature perturbation is conserved for adiabatic pertur- 
bations and thus Eq. (167|) is directly related to the amplitude of the observed primordial 
density perturbation. If the radiation and matter content in today's universe comes solely 
from the single field that dominates the final ekpyrotic phase, we will have no isocurvaure 
perturbations in an expanding universe. 

It is interesting to compare this multi-field model with a single field model that gives a 
scale- invariant spectrum of comoving curvature perturbations during collapse [HI [26J [27] . 
The single-field model requires the correct exponent for a relatively flat and positive expo- 
nential potential in order to obtain a oc |t| 2//3 , whereas the ekyprotic model only requires 
sufficiently steep, negative exponential potentials to obtain a cx \t\ 1 ^ with e ^> 1. On the 
other hand both models require fine-tuned initial conditions as it is the existence of an insta- 
bility that gives rise to the scale-invariant perturbation spectrum during collapse. In both 
models the amplitude of tensor perturbations is determined by the Hubble scale when the 
perturbations leave the horizon as tensor perturbations are then frozen on super-horizon 
scales. In the single field model the tensor metric perturbations thus acquire the same al- 
most scale-invariant spectrum [28J as the comoving curvature perturbation, with a similar 
amplitude, giving rise to a dangerously large tensor-scalar ratio which severely constrains 
the model [27]. But in the ekpyrotic scaling solution the tensor perturbations have a steep 
blue spectrum and are completely negligible at scales relevant for the cosmic microwave 
anisotropies. 

In summary, a simple ekpyrotic model with two steep, negative exponential potentials 
is capable of generating a scale-invariant spectrum of comoving curvature perturbations. A 
key ingredient is the instability of the multi-field scaling solution [T6l [TT| [121 [19] . This insta- 
bility generates a scale-invariant spectrum of isocurvature field perturbations from vacuum 
fluctuations about the scaling solution, and converts them to a scale-invariant spectrum of 
comoving curvature perturbations. We should emphasize that this conversion occurs auto- 
matically due to the dynamics of the fields in our simple model and does not require any 
change in the shape of the potential or any additional dynamics [Til [T2~[ IT3] . 

Thus ekpyrotic collapse with multiple fields can generate a scale-invariant spectrum for 
curvature perturbations in several different ways (see also [30] for a different idea). In all 
these models, we need some preceding phase that initially drives the classical background 
solution to the unstable multi-field scaling solution throughout our observable region of 
space. This problem cannot be solved within the simplest model with multiple exponential 
potentials that we considered in this paper and we would need to appeal to a more ambitious 
framework for the model such as the cyclic scenario [311 [32] to address this problem. 
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